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The Agulhas Current (AC) is the strongest western boundary current in the Southern Hemisphere and is key 
for weather and climate patterns, both regionally and globally. Its heat transfer into both the midlatitude 
South Indian Ocean and South Atlantic is of global significance. A new composite coral record (Ifaty and 
Tulear massive Pontes corals), is linked to historical AC sea surface temperature (SST) instrumental data, 
showing robust correlations. The composite coral SST data start in 1660 and comprise 200 years more than 
the AC instrumental record. Numerical modelling exhibits that this new coral derived SST record is 
representative for the wider core region of the AC. AC SSTs variabilities show distinct cooling through the 
Little Ice Age and warming during the late 18 th , 19th and 20th century, with significant decadal variability 
superimposed. Furthermore, the AC SSTs are teleconnected with the broad southern Indian and Atlantic 
Oceans, showing that the AC system is pivotal for inter-ocean heat exchange south of Africa. 



The greater Agulhas Current (AC) system near the southern tip of Africa is a key component in the global 
climate system through its role in inter-ocean heat and salt transport, thereby influencing the Atlantic 
meridional overturning circulation 1 " 3 (hereafter AMOC). Modeling experiments suggest that the increased 
export of AC waters into the Atlantic Ocean, through the so-called Agulhas leakage, results in an enhanced 
AMOC albeit with a lag of 15-30 years 4 " 7 . There is mounting evidence that the AC SSTs have increased since the 
early 1980's 8 . The recent warming is thought to be related to an increase in ocean heat transport in response to an 
increase in wind stress curl in the southern Indian Ocean trade winds 8 ' 9 . 

The AC is the strongest western boundary current in the Southern Hemisphere and releases a lot of heat to the 
atmosphere throughout the year. Relatively strong winds over the southern AC lead to strong latent heat loss from 
the surface, while further north off the east coast of South Africa, there is often convective cloud cover over its 
warm core 10 " 13 . As a result, under suitable atmospheric conditions, the AC can influence severe weather systems 
over South Africa, Botswana, Zimbabwe and Mozambique such as cut-off lows 14 , thunderstorms 15 , and mesoscale 
convective complexes 16 . 

The AC region is modulated by interannual and decadal climate modes in the Indian and Pacific Ocean, partly 
through their influence on the circulation in the South Indian Ocean 17 " 20 and also through local ocean-atmo- 
sphere interactions driven by these modes 17,21 " 23 . This strong interannual and decadal variability may be super- 
imposed on the long-term trend in AC SST that could partly be related to changes in the subtropical Indian Ocean 
gyre strength on these time scales 24-26 . The Indian Ocean gyres weaken during El Nino or positive Indian Ocean 
Dipole (hereafter IOD) events whereas they strengthen during La Nina or negative IOD events. The same holds 
true for subtropical Indian Ocean dipole (hereafter SIOD) phases 18 ' 20 . On decadal time scales, ENSO-like (El 
Nino-Southern Oscillation) decadal variability and/or the Pacific Decadal Oscillation (PDO) is thought to exert 
an influence on southwestern Indian Ocean SST and atmospheric circulation 22,27 . 

Historical gridded SST data and reconstructions reveal that the southern Mozambique Channel (MC), one of 
the source regions of the AC, has also experienced one of the highest rates in ocean warming across the southern 
Indian Ocean since 1950 28 . The instrumental SST reconstructions for the southern MC and the AC core regions 
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Figure 1 | Surface ocean connectivity between southwest Madagascar and the Agulhas Current, (a) Sea surface temperatures for February 1994 (austral 
summer) across the southwestern Indian Ocean with the major ocean currents indicated: SEC = South Equatorial Current, EMC = East Madagascar 
Current, AC = Agulhas Current, ACR = Agulhas Current return flow, MCE = Mozambique Channel Eddies. The Ifaty coral core location (yellow dot) 
and the region of dense SST observations (rectangular box) for the Agulhas Current are indicated, (b) Sea surface temperature anomalies (SSTA) relative 
to the 1961 to 1990 period for the grid box closest to the Ifaty/Tulear coral core sites (red) compared to the Agulhas Current core region data from 
HadSSTl 32 (green), ERSSTv3b 34 (blue) and HadSST3 33 (black). Note the strong co -variability of Ifaty and Agulhas core region SST on multi-decadal time 
scales. Figure la produced at http://iridl.ldeo.columbia.edu/ and used with permission (data from ref. 44). 



show considerable multidecadal variability (Fig. 1). A single 334 year 
long coral geochemical proxy record from Ifaty reef off southwestern 
Madagascar in the southern MC also revealed an increase in SST after 
the 1970's 20 . This coral record revealed a non- stationary relationship 
of local southern MC SST with ENSO, and a link with the Pacific 
Decadal Oscillation 27 . These single core observations, however, are 
limited and may not represent the wider southern MC or the AC core 
region as yet. 

The aims of this paper are to produce a new composite coral record 
to obtain a SST record representing SST changes in the wider south- 
ern MC and the AC core regions. Here, we present a new three-core 
coral composite oxygen isotope record from the Ifaty and Tulear 
coral reefs off SW Madagascar, which covers SST variability for the 
past 334 years (1660-1994). Through modeling, we investigated 
whether the southern MC SST data are representative for the down- 
stream AC core region. Finally, we aim to infer cross -ocean relation- 
ships in the Southern Hemisphere using SST records, elucidating the 
pivotal role of the AC region at the cross roads between the Indian 
and Atlantic Oceans. 



Results 

Connectivity between the Ifaty coral site and the Agulhas Current 
core region, at 32 °S. First, we tested connectivity in ocean 
properties, namely SST and salinity (SSS), between the Ifaty coral 
site and the Agulhas core region at 32°S, the so-called CAP 
(Lagrangian virtual floats capture area) area in Figure 2. We use 
the 5-daily output from a 1948 to 2007 hindcast of a 1/4°, eddy- 
permitting, basin-scale ocean model 29 (Fig. 2). To establish 
connectivity between the Ifaty site and the downstream AC core 
region, we correlate the detrended monthly-mean SSS and SST 
time series from Ifaty with those extracted across the AC source 
regions and western Indian Ocean boundary (Fig. 2a, b). Our 
results suggest that SSS and SST across most of the regions of the 
southern MC and southwestern Indian Ocean boundary are 
correlated. Further, we find a strong SST and SSS relationship 
between the Ifaty site and AC core region, shown by the same 
inter-annual surface variability with zero lag (Fig. 2c). 

To assess the time scales of horizontal ocean transport (advection) 
between the Ifaty site and the AC at 32°S, Lagrangian virtual floats 
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Figure 2 | Connectivity of sea surface temperature (SST), salinity (SSS) and near surface Lagrangian derived advection. (a) Regional correlation with 
the detrended monthly-mean IFATY surface salinity (SSS) time-series for 1948-2007 (shading, p < 0.05). Correlations are plotted at maximum lag, 
with lag values shown as contours, (b) as in a) but for SST. (c) lag correlation for SST and SSS between Ifaty and the Agulhas Current at 32°S 
(represented by the IFA point and AC32 box in panel (a), (d) Advective connection between Ifaty and the Agulhas Current, as derived from Lagrangian 
floats released at RE (panel b), and collected at CAP (panel b). Figures produced in MATLAB. 



are deployed in the surface waters (z < 100 m) at the former location. 
Figure 2d shows that 76.3% of the particles released reach the AC in 3 
years, with a median transit time between the two sites of 240 days. 
The advective time scale is quite long, whereas the lag between the 
two SST and SSS signals is near zero (Fig. 2c). This suggests two 
possibilities: 1) that both sites happen to 'receive' the same SST 
and SSS signal through ocean transport from upstream, or 2) that 
the signals are driven by large-scale atmospheric processes, influ- 
encing both the wider MC area and the AC core region. 
Nevertheless, our model data confirm the physical link in SST 
between the MC and the AC core regions and return flow region. 
The modeling experiments show that the spatial and temporal cor- 
relation of SSTs and SSS of the wider Agulhas region is high, includ- 
ing the source areas between the coral sites in the southern MC and 
the AC core region, for instance in the CAP area at 32°S. 

Coral-based SST reconstruction for the southern Mozambique 
Channel and the Agulhas Current core region. The new SST 

reconstruction is based on three Porites coral oxygen isotope 
(5 ls O) records at annual resolution from Ifaty and Tulear coral 
reefs off southwestern Madagascar, 43°E, 23°S, covering the past 
334 years (Fig. 3a; Fig. SI). The new coral composite comprises 
three corals from 1905 to 1994, two from 1881 onwards and a 
single coral covering 1660 to 1880 (Fig. 3b). The coral proxy data 
and measurement procedures are described in detail in the Methods. 
In brief, all proxy records were centered by removing the 1961-1990 
mean. A composite annually resolved coral temperature record was 



then constructed by (1) converting each 5 ls O record to temperature 
units, (2) calculating the arithmetic mean of the coral SST records 
from each site, and (3) averaging the mean SST records from both 
sites (Fig. SI). For the proxy- temperature conversion, we use the 
mean of the published 5 ls O- temperature relationship 30 ( — 0.2%o/°C 
for 5 ls O). This has the advantage that our coral temperature 
reconstruction does not depend on linear regression with instru- 
mental SST data, which has rather large errors 31 . The composite 
chronology extends from 1660 to 1994 (Fig. 3). 

The composite SST record indicates strong multidecadal variabil- 
ity throughout the 334 years. The composite mirrors the interannual 
and decadal variability of AC core region SST since 1854 in ERSST, 
since 1850 in HadSST3 and since 1870 in HadSSTl (Figure 3). The 
coolest period in SW Madagascar and the AC core region on record 
pre- 1900 is observed between 1675 and 1720 and includes the Late 
Maunder Minimum (1675-1710), the coolest period during the Little 
Ice Age. The cool period shows strong interannual variability. The 
warmest periods on record pre- 1900 for both SW Madagascar and 
the AC core region are observed between 1660-1670, 1770 to 1805, 
and 1870 to 1900. The warming rate between 1720 and 1800 was 
rapid at 0.13 ± 0.02°C/decade. The warm period between 1870 to 
1900 is also very prominent for both SW Madagascar and the AC 
core region. This period is punctuated by strong interannual warm 
anomalies with cooler years in between. Similarly strong interannual 
variability is also observed for the 1770 to 1805 warm period. SST was 
lower in the early 20th century and increased markedly after 1970 for 
both SW Madagascar and the AC core region. 
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Figure 3 | Sea surface temperature reconstruction for the Agulhas Current region, (a) Coral composite SST anomaly reconstruction for southwestern 
Madagascar (grey shading shows 2 standard error) compared to Agulhas Current core region data from HadSSTl 32 (green), ERSSTv3b 34 (blue) and 
HadSST3 33 (black), (b) Time coverage of individual coral core records and (c) Same as (a), yet for filtered data (decadal) with a cutoff of 8 years. All SST 
anomalies computed relative to the 1961 to 1990 period. 



Validation of SST trends in the instrumental and coral SST records. 

Absolute warming rates for our coral index and instrumental SST have 
been estimated using ordinary least squares (OLS) regression, which 
mainly depends on the 5 ls O-SST conversion. We have estimated the 
warming rates obtained using published 5 ls O-SST relationships 30 with 
the mean, upper and lower estimates as uncertainty bounds for the 
regression (—0.2 ± 0.02%o/°C). We use the period since 1900 because 
SST products are less reliable prior to that date. 

We found a steady warming in the coral composite SST record 
since 1900 (1.07 ± 0.1°C), with most of the increase in the warming 
rates after 1970 (1.03 ± 0.3° C). The magnitude of the warming 
inferred from the corals is 0.12 ± 0.02°C/decade since 1900 and 
ranges between 0.13°C and 0.11°C/decade (using the upper and 
lower limits of 5 ls O-SST relationships). We speculate that much of 
the observed spread in proxy- SST relationships reflects statistical 
uncertainties since the error of the OLS regression obtained using 
the mean slopes equals the spread of the warming rates estimates 
using the upper and lower limits of proxy-SST relationships. The 
coral composite SST reconstruction has a higher correlation with 
instrumental SST for both the Ifaty-Tulear and Agulhas Current core 
regions than any single coral record (Table 1; Table SI and S2). The 



warming trend after the mid-1970's in the coral composite record 
agrees substantially better with instrumental SST than the previously 
published Ifaty-4 coral record alone 20 . 

To validate the trend in our coral composite SST record, we 
extracted SST records for the Ifaty-Tulear region (centred at 43 °E, 
23°S) and the AC core region (26°-28°E, 34°-36°S) from three grid- 
ded SST products for the annual average computed from March to 
February: HadSSTl 32 from 1870, HadSST3 33 from 1850, and 
ERSSTv3b (hereafter ERSST 34 ) from 1854 (Fig. 3). All SST products 
are based on ICO ADS data 35 , but were constructed with different 
strategies for historical bias correction and gridding. SST data for the 
Ifaty-Tulear region is based on extremely sparse observations in the 
ICO ADS 35 data (Fig. S2). The new 5° X 5° gridded HadSST3 data set 
is based on a comprehensive reanalysis of the data and metadata in 
the ICOADS database and is bias-corrected with improved error 
estimates 33 . We identified the gridbox in HadSST3 along the 
Agulhas Current (24°-29°E, 31°-36°S) which is the most complete 
back to 1850 for the greater Agulhas Current region (Fig. 1). 

The long-term warming rates indicated by the SST products for 
the Ifaty-Tulear region since 1900 clearly differ, yet agree on the 
strong post 1970 warming. The warming rates per decade since 
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Table 1 | Correlations of detrended, mean annual (March to Feburary) reconstructed SST (IF comp = Ifaty/Tulear coral composite SST) with 
instrumental SST at Ifaty (IF) and for the Agulhas Current (AC) core region. Hadl = HadSSTl 32 , Had3 = HadSST3 33 , ERSST 34 . All 
correlations significant above the 1% level, despite IF comp with IF Had3 (5%). Correlations computed for maximum number of years in 
each dataset taking into account the degrees of freedom for each correlation pair. Ifaty HadSST3 suffers from sparse observation pre-1 920, 
while AC HadSST3 is complete to 1 850. Correlations computed at http://climexp.knmi.nl/ 42 
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1900 are 0.054 ± 0.008°C (ERSST 34 ), 0.081 ± 0.01 (HadSSTl 32 ) and 
0.043 ± 0.01 (HadSST3 33 since 1909). The absolute warming 
between 1970 to 2005 ranges between 0.39 ± 0.15°C (HadSST3), 
0.42 ± 0.1°C (ERSST) and 0.71 ± 0.12°C (HadSSTl). 

Our coral-based SST reconstruction tracks instrumental AC core 
region SST since 1854 in ERSST, since 1870 in HadSSTl and since 
1850 in HadSST3 (Fig. 3; Tab. 1). We found significant correlations 
for detrended coral composite SST with AC core region HadSSTl 
(r = 0.57, P < 0.001, DF = 124), HadSST3 (r = 0.59, p < 0.001, 
DF = 142) and ERSST (r = 0.54, p < 0.01, DF = 140) for the 
maximum number of years in each SST dataset (Tab. 1). The correla- 
tions were still significant for an early verification period between 1850 
to 1910 with AC HadSST3 (r = 0.62, p = 0.0005, DF = 59), although 
lower with HadSSTl (r = 0.54, p = 0.005) and ERSST (r = 0.4, 
p = 0.03). The correlation between the coral composite and AC core 
region SST is higher than with Ifaty- Tulear region SST from all SST 
products (Tab. 1; Tab. SI). This is most probably due to the sparse to 
non-existent observations for the Ifaty-Tulear region (Fig. S2). 

All records, including the coral composite SST, show multidecadal 
to centennial variability in SST since at least 1870 with two warm 
periods between 1870 to -1900 and 1970 to 2005, bracketed by 
cooler SST in between (Fig. 3). Interannual warm peaks are observed 
during the 20th century, the most prominent between 1939-1942 
and the late 1950's to early 1960's. 

Relationships with regional and large-scale temperatures. Our 

new SST reconstructions also show excellent correlation with an 
annual mean South African air temperature reconstruction (AT) 



from a stalagmite record 36 over the entire 334 year period (Fig. 4). 
The correlation is strongest on decadal time scales applying a Loess 
filter with a cutoff at 8 years (r = 0.59, p = 0.008, N = 308; Fig. 3c). 
Both, the coral composite SST and stalagmite AT indicate the period 
between 1690 to 1740 as the coolest on record. In addition, both 
records indicate relatively high temperatures during the late 18th 
and 19th century, a cool period in the early 20th century followed 
by a warming towards the end of the 20th century. 

To assess the correlations of the coral composite with global SST 
on decadal time scales, we applied a Loess low-pass filter with a cutoff 
at 8 years. Spatial correlations between coral composite SST and 
global SST, from HadSSTl, revealed statistically significant correla- 
tions (p < 0.05) on decadal time scales (Fig. 5) The highest correla- 
tions are found in the AC region and across the southern Indian 
Ocean between 20 to 60°S stretching towards the western and south- 
ern Australian coast (Fig. 5). Positive correlations also emerge off the 
eastern coast of South America and the (sub)tropical Atlantic. Other 
teleconnected regions with positive correlations are the northern 
Indian Ocean, the western and northeastern Pacific. 

We also performed an EOF (Empirical Orthogonal Function) 
analysis of detrended mean annual HadSSTl in the southern 
Indian Ocean (10-40°S, 20-130°E) to reveal the dominant SST pat- 
terns. EOF1 of HadSSTl showed positive loading across the southern 
Indian Ocean between South Africa and western Australia. The PCI 
time series of this EOF1 pattern resembles the multidecadal SST 
record for the Agulhas Current region and the coral composite 
SST (Fig. S3). Our new coral composite SST is strongly related with 
PCI showing significant correlations (p < 0.001; Suppl. Tables S4 to 
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Figure 4 | Reconstructions for sea surface and air temperature for the Agulhas Current and South Africa. Coral composite SST anomaly reconstruction 
for southwestern Madagascar (grey shading shows 2 standard error) compared to an air temperature reconstruction from stalagmite T7 from Cold 
Air Cave (24°S, 29°E, 1420 m above sea level) in South Africa 36 (blue shading shows 2 standard error). Note the gap in the stalagmite record between 1954 
and 1973. Both SST and air temperature anomalies computed relative to the 1961 to 1990 period. 
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Figure 5 | Global teleconnections of Agulhas Current region sea surface temperature, (a) Spatial correlation of coral-derived SST reconstruction with 
global HadSSTl 32 between 1870 to 1994 at decadal time scales (>8 years) and (b) same as (a), yet for Agulhas Current SST from HadSSTl 32 . Correlation 
with p < 0.05 colored. Figure produced at http://climexp.knmi.nl/ 42 and used with permission. 



S6 and Fig. S3). EOF2 shows a dipole pattern in SST between the 
southwestern Indian Ocean east of Madagascar and the southeastern 
Indian Ocean offshore of Indonesia/NW Australia (Fig. S3). This 
pattern is related to ENSO. PC2 is not correlated with the coral 
composite when the entire period from 1870 to 1994 is considered. 
However, after 1950 the correlation is significant (r = 0.44, p < 0.01, 
DF = 40). This finding is consitent with Zinke et al. (2004) 20 who 
showed evidence for a non- stationary relationship of southern MC 
SST with ENSO. EOF3 showed a dipole pattern in SST between the 
tropical and subtropical/mid-latitude Indian Ocean. EOF3 com- 
prises the well-known subtropical dipole modes that can be found 
across all southern ocean basins, especially in austral summer. The 
PC3 time series of this EOF3 pattern is positively correlated with our 
coral composite SST (Suppl. Table S7and Fig. S3). This finding is 
consistent with Zinke et al. (2004) 20 , who showed evidence for a 
correlation of southern MC SST with the subtropical dipole mode 
of the southern Indian Ocean. 

Discussion 

Our new annually resolved coral composite SST reconstruction 
tracks instrumental AC core region SST since the late 19th century 
observed in the three SST products: ERSST, HadSSTl and HadSST3. 
Based on this strong relationship, we are able to resolve historical SST 
variability and long-term changes in the southern MC and the AC 
core region back to 1660. We observed strong multidecadal to cent- 
ennial-scale oscillations in SST over the past 334 years in the coral 
composite SST and since at least 1850 in three SST datasets. We 
confirmed the period between 1670 and 1720 as the coolest period 
on record, which suggests that transport of warm water from the 
tropical Indian Ocean towards the AC was diminished at this time. 
We find warm excursions between 1770 and 1800, the late 19th 
century and post 1976 where warm water transport to the AC was 
at its peak. The current warming after 1976 to 1994 is not unpreced- 
ented since similar warming rates were observed from 1720 to 1800 
following peak cooling during the Late Maunder Minimum during 
the Little Ice Age. Yet, instrumental SST beyond 1994 show that the 
warming of the past 20 years was exceptional in the context of the 
past 334 years 8 (Figs. 1 and 3). However, changes in 5 ls O of seawater 
through ocean advection and/or precipitation- evaporation changes 
could have contributed to lower coral 5 ls O and therefore a warm bias 
during the late 18th century. A previous study coupled (bi) monthly 
resolved Ifaty 5 ls O and Sr/Ca SST proxies for the a single core (Ifaty- 
4) and suggested that variations in 5 ls O of seawater might have 
influenced the SST estimates during the late 18th and 19th centur- 
ies 20 . These periods stand out as particularly warm decades in our 
new reconstruction. Since our new reconstruction pre- 1882 is still 



based on the Ifaty-4 record only, a 5 ls O seawater contribution cannot 
be excluded. Nevertheless, the three instrumental SST datasets agree 
with the coral composite SST in showing a warm period with strong 
interannual oscillations pre- 1900. 

In addition, the recently published air temperature reconstruction 
for southern Africa from a stalagmite record 36 agrees with the coral 
composite in showing higher temperatures in the region during the 
late 18th and 19th century (Fig. 4). A rainfall reconstruction for 
southern Africa 37 revealed that the 19th century was wetter than 
the 20th century which could be partly related to relatively warm 
SST in the AC region during the early and late 19th century observed 
in our SST reconstruction. Using observations and AGCM experi- 
ments, Reason and Mulenga (1999) 38 showed that SSTs in the AC 
region are related to South African rainfall on interannual and multi- 
decadal time scales. This adds confidence to the reality of a relatively 
warm late 19th century in the AC region, only matched and/or 
exceeded in the period after 1976 (Figs. 1, 3 and 4). 

Our reconstruction is in concert with typical spatial teleconnec- 
tion pattern of the AC region SST with global SST on interannual and 
decadal time scales. We did not find an ENSO spatial correlation 
pattern, which agreed with the weak or non- stationary ENSO cor- 
relation on seasonal and mean annual time scales with instrumental 
SST for the AC region 20 . Yet, we found significant correlations with 
the western and northeastern Pacific SST. This pattern resembles to 
some extent the ENSO -like decadal or Pacific Decadal Oscillation 
pattern in SST 27 . The strongest relationship was found with the 
southern Indian Ocean and off the east coast of South America. 
The strong correlations with the southern Indian Ocean are found 
along the AC return flow. The teleconnections with the South 
Atlantic are found along the advective pathways of the Agulhas 
leakage 7,39 . 

The strongest SST teleconnections across the southern Indian and 
Atlantic Ocean point to a dominant role of the AC system for SST 
variability across the subtropical to mid-latitude Indian and Atlantic 
Ocean (Fig. 5). A recent study of mid-latitude islands across the 
southern part of all ocean basins also revealed that the Agulhas 
system is likely a dominant driver of climate variability for the 
mid- to high latitude Indian Ocean islands 40 . The recent spatial 
warming trend pattern in the southern Indian Ocean is characterised 
by two warming centers, one in the southwestern and the other in the 
southeastern Indian Ocean (Fig. S3). This pattern strongly resembles 
EOF1 of HadSSTl for the southern Indian Ocean (Fig. S3). PCI of 
this EOF1 showed multidecadal SST variability and no linear warm- 
ing trend (Fig. S3). Thus, the two centers of warming most possibly 
underwent similar multidecadal SST changes and the recent pattern 
is part of this natural multidecadal variability. The mechanism 
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driving this multidecadal variability might be related to large-scale 
wind forcing across the southern Indian and Atlantic Ocean 17 ' 41 . Lee 
et al. (201 1) 39 showed that the increased heat content across the 
southern Atlantic since the 1950's was largely related to an increase 
in inter-ocean heat exchange with the Indian Ocean. The warming in 
the South Atlantic after 1950 was probably reinforced by a warm 
Agulhas leakage and its westward and northward transport into the 
Atlantic. Similar warming of the South Atlantic should have 
occurred during the late 18th and 19th century according to our 
AC core region SST reconstruction. This assumption needs to be 
tested with coupled atmosphere-ocean models. 

In conclusion, our new coral composite SST record for the 
Agulhas Current system showed strong multidecadal SST variability 
in this globally important ocean current and that the warming over 
the last three decades is not unprecedented in the context of the 
multi- centennial record. 

Methods 

Coral core collection and sampling. Coral cores from massive Pontes sp. at Ifaty and 
Tulear reefs were collected in October 1995 during the European Union TESTREFF 
program from the Ifaty- Ranobe lagoon and the Great Barrier of Tulear (southwest 
Madagascar 20 ). The Ifaty and Tulear coral reef sites are described in detail in Zinke 
et al., 2004 20 and Bruggeman et al., 2012 43 , respectively. Core Ifaty-4 (4.06 m length), 
core Ifaty- 1 ( 1 .93 m length) and Tulear-3 ( 1 .80 m length) were obtained from a depth 
of 1.1 m, 1.8 m and 0.6 m below mean tide level. The average growth rate of core 
Ifaty-4 was 0.99 ±0.15 cm per year, whereas Ifaty- 1 averaged 1.28 ± 0.24 cm and 
Tulear-3 1.54 ± 0.25 cm per year. 

All cores were sectioned to a thickness of 7 mm and slabs were cleaned in 10% 
hydrogen peroxide for 48 h to remove organic matter at GEOMAR Kiel. Then, slabs 
were rinsed several times with demineralized water and dried with compressed air. For 
complete removal of any moisture within the coral skeleton the sample was put into 
an oven for 24 h at 40°C. The slabs were X-rayed to reveal annual density banding. 

A high resolution profile for stable isotope analysis on core Ifaty-4 was drilled using 
a computer-controlled drilling device along the growth axis as observed in X-ray- 
radiograph-positive prints 20 . Subsamples were drilled at a distance of 1 mm for the 
years 1995-1920 and 2 mm for the older part of the core; the drilling depth was 3 mm 
using a 0.5 mm dental drill at 1000 rpm. The 1 or 2 mm sample spacing provides 
approximately monthly or bimonthly resolution, respectively. Cores Ifaty- 1 and 
Tular-3 were sampled at annual resolution along the major growth axis following the 
density pattern from summer to summer in any given year established from X-ray- 
radiograph-positive prints. 

Analytical procedures 5 18 0. The high-resolution samples of core Ifaty-4 were 
reacted with 100% H 3 P0 4 at 75°C in an automated carbonate reaction device (Kiel 
Device) connected to a Finnigan MAT 252 mass spectrometer (University Erlangen). 
Average precision based on duplicate sample analysis and on multiple analysis of NBS 
19 is ±0.07%o for 5 ls O (la). The annual samples were reacted with 100% H 3 P0 4 at 
75°C in an automated carbonate reaction device (Kiel Device) connected to a 
Finnigan MAT 252 mass spectrometer at the VU University of Amsterdam. Average 
precision based on duplicate sample analysis and on multiple analysis of NBS 19 is 
±0.08%ofor5 18 O(la). 

Proxy data treatment. We used the already published, (bi) monthly resolved Ifaty-4 
coral 5 ls O time series from 1660 to 1994 20 . The high resolution of the Ifaty-4 coral 
5 18 0 enabled us to compute a precise annual chronology averaged between March to 
February. We used the Ifaty-4 core as our best dated reference time series to ensure 
that the yearly sampled chronologies of Ifaty- 1 and Tulear-3 aligned well. The three 
proxy records were first centered by removing the 1961-1990 mean. A composite 
annually resolved coral temperature record was then constructed by (1) converting 
each proxy record to temperature units, (2) calculating the arithmetic mean of the 
coral SST records from each site, and (3) averaging the mean records from both sites. 
For the proxy-temperature conversion we use the mean of the published § 18 0- 
temperature relationship ( — 0.2 ± 0.02%o/°C for 5 18 0 30 ). This resulted in a time series 
of relative SST changes at annual resolution against the 1961 to 1990 mean. 

We estimated the uncertainties of the coral composite SST reconstruction fol- 
lowing the method of Nurhati et al. (201 1) 31 for relative SST reconstructions. The 
errors are displayed in Figure 3 as grey shaded envelopes for each individual year. 
Uncertainty in the relative annual mean SST reconstruction takes into account errors 
associated with (i) the analytical precision of the 5 ls O measurements (±0.4°C), and 

(ii) the calibration slope of the 5 ls O-SST calibration (±0.1°C). Taken together, the 
annual mean coral 5 ls O-derived SST error is ±0.41°C (la), quadratically combining 
terms (i)-(ii). 

Uncertainty estimates in the relative SST reconstruction trend was estimated 
taking into account errors associated with (i) the slope error of the SST trend 1900 to 
1994 (±0.16°C), (ii) the analytical precision of the § ls O measurements (±0.4°C), and 

(iii) the calibration slope error of the 5 ls O-SST calibration (±0.1°C). Taken together, 
the twentieth-century coral 5 ls O-derived SST trend error is 0.44°C, quadratically 
combining terms (i)-(iii). 
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